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Abstract 


Fosberg,  Michael  A.,  William  E.  Marlatt,  and  Lawrence  Krupnak. 
1976.  Estimating  airflow  patterns  over  complex  terrain.  USDA 
For.  Serv.  Res.  Pap.  RM-162,  16  p.  Rocky  Mt.  For.  and  Range 
Exp.  Stn.,  Fort  Collins,  Colo.  80521. 

A  simple  one-layer  model  of  atmospheric  boundary  layer  flow 
was  developed  for  use  in  complex  terrain.  The  model  was  derived 
through  simplification  of  the  fundamental  Navier-Stokes  flow  equa- 
tions. This  simplification  was  gained  by  neglecting  advection  terms 
in  the  equations  of  motion  and  by  assuming  an  impulse  solution. 

Although  the  equations  do  not  describe  all  flow  characteristics, 
the  resultant  solution  describes  a  diagnostic  model  of  the  vector 
flow  field.  It  requires  much  less  data  than  traditional  approaches, 
and  therefore  can  be  used  as  an  estimator  of  wind  patterns  in  areas 
where  dense  observational  networks  are  not  economically  feasible. 

The  intended  primary  uses  for  this  model  are  in  providing  wind 
fields  for  fire  behavior  prediction  and  in  evaluation  of  pollution  trans- 
port patterns. 
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Estimating  Airflow  Patterns  Over 
Complex  Terrain 


Michael  A.  Fosberg,  William  E.  Marlatt, 
and  Lawrence  Krupnak 


Introduction 

Land  managers  are  increasingly  faced  with 
problems  requiring  quantitative  estimates  of  air- 
flow patterns  over  complex  terrain.  Use  of  the 
wind  information  may  range  from  potential  air 
quality  analysis  for  land  use  planning  to  predict- 
ing spread  of  actual  wildfires.  Areas  in  which  such 
analysis  is  required  are  typically  mountainous. 
Nearly  always,  evaluation  is  severely  hampered 
by  a  single  problem— lack  of  data,  particularly 
wind  data. 

For  these  reasons,  a  method  of  determining 
boundary  layer  flow  patterns  is  needed  as  an  alter- 
native to  expensive  data  networks. 

Boundary  layer  flow  over  complex  terrain 
cannot  currently  be  numerically  modeled  with 
complete  fluid  and  thermodynamic  equations. 
Accuracy  and  quantity  of  the  initial  data  fields 
required  for  complete  solution  only  shift  the  prob- 
lem back  to  need  for  expensive  data  networks 
(Fosberg  1967,  1969).  While  these  costs  and  efforts 
may  be  worthwhile  in  special  locations  and  for 
specific  problems,  flexibility  in  geographic  loca- 
tion and  timeliness  of  solutions  are  sacrificed.  We 
feel  that  an  alternative  procedure  is  required  for 
wind  pattern  analysis  which  provides  geographic 
flexibihty  and  does  not  require  sophisticated, 
expensive  data  networks.  Such  solutions  cannot 
include  all  atmospheric  processes,  however.  The 
governing  equations  must  be  simplified  and  some 
important  effects  must  be  represented  by  param- 
eterized analogs. 

A  model  of  this  type  is  based  on  tradeoffs 
between  simplified  solutions  and  preservation  of 
realistic  results.  Several  models  of  thermal  and 
terrain-induced  perturbation  of  near-ground  flow 
have  been  developed  recently.  Anderson  (1971) 


and  Lantz"*  developed  single  layer  models  of  two- 
dimensional  potential  flow  over  terrain  features. 
Anderson  (1971),  Fosberg  et  al.  (1972),  and  Fos- 
berg^ examined  the  thermally  driven  flow  dis- 
turbances. These  studies  were  all  based  on  the 
philosophical  argument  that  major  simplification 
of  the  equations  was  acceptable,  and  that  the  dis- 
turbed flow  field  could  be  superimposed  on  a  mean 
flow  field.  While  superficially  attractive,  these 
preliminary  models  did  not  evaluate  the  conse- 
quences of  the  simplifications,  nor  did  they  provide 
quantitative  arguments  supporting  neglect  of 
specific  processes. 

Each  of  the  models  describes  a  single  layer  of 
finite  depth,  depicts  a  quasi-horizontal  slab,  and 
treats  body  forces  as  an  impulse  acceleration.  The 
models  also  require  the  flow  to  be  statically  stable, 
either  as  an  incompressible  fluid  or  through  the 
hydrostatic  assumption.  Superposition  of  dis- 
turbed flows  on  the  mean  flow  also  implies  Umita- 
tions  on  the  processes  that  can  be  modeled.  Im- 
pulse accelerations  tacitly  imply  that  a  model  is 
intended  as  a  diagnostic  tool  rather  than  a  time- 
dependent  predictive  model.  Also,  these  models 
generally  neglect  or  highly  linearize  the  advection 
terms  in  the  dynamic  equations. 

The  highly  idealized  treatment  of  the  advec- 
tion term  is  probably  the  most  serious  of  the 
assumptions.  Justification  is  provided  partially 
by  examination  of  the  energy  integrals  in  Fos- 

*Lantz.  R.  B.  1972.  Application  of  a  tliree  dimensional  numer- 
ical model  to  air  pollution  calculations.  Preprint  of  paper  presented 
at  65th  Annu.  I^eet.,  Air  Pollution  Control  Assoc,  Miami  Beacfi, 
Fla..June  18-22,  22  p. 

^Fosberg,  Michael  A.  (in press).  New  technology  for  deter- 
mining atmospheric  influence  on  smoke  concentrations. 
Proc.  Int.  Symp.  on  Air  Quality  and  Smoke  from  Urban  and 
For.  Fires.  1973.  lip. 
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berg's  (1967,  1969)  analysis  of  flow  over  single 
ridges.  Conversion  of  potential  energy  to  kinetic 
energy  dominates  the  kinetic  energy  budget.  The 
advective  redistribution  of  kinetic  energy  and  the 
viscous  dissipation  are  about  an  order  of  magni- 
tude smaller.  While  these  energy  budget  compari- 
sons do  not  fully  justify  neglect  of  the  advective 
processes,  they  do  show  that  the  assumption  does 
not  invaUdate  the  models,  but  only  restricts  inter- 
pretation of  the  results. 

In  general,  this  class  of  airflow  models  ex- 
cludes the  small-scale  physical  situations  jn  which 
transient  or  highly  interactive  processes  such  as 
gravity  waves  or  unstable  thermal  convection 
dominate  the  flow  field.  Also,  because  of  the 
restrictions  imposed  on  the  advection  terms, 
downwind  wake  effects  are  not  properly  repre- 
sented. 

Use  of  these  simplified  models— rather  than 
complex  and  complete  solutions  or  as  a  substitute 
for  extensive  field  studies— becomes  attractive 
when  extensive  data  are  not  available,  or  in  land 
use  planning  where  alternatives  may  be  great,  the 
time  span  for  planning  short,  and  the  accuracy  or 
costs  of  more  physically  precise  solutions  and 
studies  unwarranted. 

The  model  described  in  this  paper  is  intended 
to  meet  these  applied  objectives.  This  model  is 
based  on  simplified  solutions  of  the  complete 
equations.  It  is  philosophically  identical  to  the 
preliminary  models,  but  combines  more  of  the 
physical  processes. 


Structure  of  the  Model 

The  mathematical  model  is  based  on  the 
terrain-induced  flow  processes  derived  by  Ander- 
son (1971)  and  the  thermally  induced  flows  de- 
scribed by  Fosberg  et  al.  (1972).  In  addition,  the 
frictional  terms  are  included  in  the  model.  Ander- 
son's model  was  based  only  on  a  form  of  the  diver- 
gence equation,  while  the  Fosberg  et  al.  (1972) 
model  was  based  only  on  the  vorticity  equation. 
The  model  developed  here  is  developed  from  both 
equations. 

The  advection  terms  were  neglected  in  both 
previous  equations.  The  consequences  of  this 
assumption  are  that  the  equations  were  simplified 
at  the  expense  of  losing  a  complete  physical  de- 
scription of  the  dynamic  processes.  Thermally, 
frictionally  and  terrain-induced  flows  were  super- 
imposed on  a  background  potential  flow  across 
the  computational  boundaries.  Each  of  these 
disturbances  was  assumed  to  take  place  only 
within  the  computational  area.  Superposition  of 
disturbances  on  the  background  flow,  after  ne- 


glecting advection  terms,  is  physically  analagous 
to  an  impulse  acceleration.  The  diagnostic  nature 
of  the  model  was  achieved  by  allowing  the  impulse 
to  act  over  a  small,  finite  time  interval  defined  by 
phase  velocities  of  disturbances.  Flow  following 
the  ground  surface  was  obtained  by  a  coordinate 
transform.  A  rigid  upper  surface  above  the  ter- 
rain was  assumed  in  order  to  define  the  top  of  the 
atmospheric  slab. 

Procedures  for  solving  the  equations  involve 
serial  approximations  which  superimpose  a  new 
physical  effect  on  the  previous  solution.  The  first 
step  in  the  solution  was  to  transfer  the  large-scale 
background  wind  into  a  terrain-induced  modifica- 
tion of  the  throughflow.  This  step  provided  a  local 
throughflow  wind. 

Next,  thermal  and  frictional  modification  of 
the  vorticity  and  divergence  were  introduced. 
These  changes  were  superimposed  on  the  terrain- 
induced  flow.  Finally,  stream  functions  and  veloc- 
ity potentials  were  calculated  so  that  the  wind- 
speed  and  direction  could  be  defined  at  all  interior 
points. 

Throughflow  Representation 

Total  airflow  through  the  computation  field 
was  required  to  satisfy  conservation  of  mass. 
Since  the  lower  surface  of  the  layer  was  mountain- 
ous, provisions  had  to  be  made  for  local  channeling 
and  deflection.  Throughflow  characteristics  were 
calculated  from  the  steady-state  mass  continuity 
equation  as  outlined  by  Anderson  (1971).  Vertical 
integration  and  use  of  the  anelastic  approxima- 
tion permitted  local  velocities  to  be  approximated 
from  the  overall  wind  field  and  terrain  features. 
The  steady-state  continuity  equation  is:^ 


In  this  equation,  V  is  the  horizontal  gradient 
operator,  p  is  density,  v  is  the  velocity  vector, 
and  w  is  the  vertical  motion,  positive  upward. 

Vector  coordinates  used  in  the  mathematics 
are  west  to  east  in  the  x  direction,  south  to  north 
in  the  y  direction,  and  vertically  as  z. 

The  anelastic  approximation,  where  kinematic 
influences  on  density  are  negligible  compared  to 
thermodynamic  influences,  was  first  applied  to 
the  equation: 

-  8(pw) 

^See  appendix  for  a  complete  description  of  notations  used 
in  the  equations. 
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Bars  over  a  variable  indicate  a  mean  value  for  the 
slab.  The  equation  was  then  integrated  over  the 
depth  of  the  layer  between  the  ground  (h)  and  the 
top  of  the  slab  (H). 


3(pw) 
8z 


/•H  /•  H 

pi     V»v  dz  =  -  / 

Horizontal  divergence  was  defined  as  d  =  V  •  y 
so  that  the  integral  equation  was  solved  as: 

p"  6"  (H-h)  =  -pw(H)  +  pw(h) 

through  the  mean  value  theorem.  A  rigid  lid  was 
assumed,  so  pw(H)  became  zero.  Vertical  motion 
near  the  ground  was  assumed  to  be  dominated  by 
orographic  lifting.  The  horizontal  divergence  of 
the  mean  wind  field  then  became: 


Vh 


(1) 


-~  (H-h) 


Terrain  modification  of  the  large-scale  flow  could 
then  be  calculated  from  the  channeling  and  deflec- 
tion of  this  flow,  and  could  provide  a  first  approxi- 
mation of  the  local  wind. 


Thermal  and  Frictional  Influences 

Modification  of  the  wind  field  produced  by 
thermal  acceleration  and  retardation  required 
invoking  the  assumptions  associated  with  the 
advection  terms  of  the  complete  equations.  Fric- 
tional influences  were  included  as  spatial  varia- 
tions of  idealized  and  parameterized  solutions 
relating  surface  characteristics  to  the  wind  profile. 
The  fundamental  procedure  was  to:  (1)  derive  the 
complete  divergence  and  vorticity  equations, 
(2)  judiciously  neglect  the  advective  terms,  (3) 
allow  full  influence  of  the  thermal  terms,  and 
(4)  utilize  classic  frictional  solutions. 

The  procedure  required  for  inclusion  of  these 
body  forces  was  to  define  the  thermal  and  fric- 
tional solutions  in  terms  of  the  Navier-Stokes 
equations.  The  first  step  was  to  express  the  diver- 
gence equation  utilizing  these  assumptions.  Diver- 
gence was  obtained  through: 

^'dF  "  It  ^  V.(-Cp0Vp^)  -f  V-(KV%) 

where  9  is  the  potential  temperature,  p^  is  the 
pressure  variable  defined  as  p^  =  TId,  and  K  is 
the  turbulent  viscosity.  The  friction  term  was 
simplified  to: 


\8z   /  \8x8z        9y8z  / 


since  the  horizontal  dissipation  terms  are  small 
compared  to  the  vertical  terms  when  expressed  in 
finite-difference  mesoscale  models.  A  parameter- 
ized form  of  the  dissipation  was  obtained  from 
the  logarithmic  wind  profile  near  the  ground.  The 
wind  profile  for  an  adiabatic  surface  layer  (Sutton 
1953)  is: 


3u. 
9z 


— r-  and  — 


(2) 


kZ 


9z  kZ 


where  u^  and  v^  are  the  friction  velocities  and  k  is 
the  von  Karman  constant.  These  two  functions 
were  differentiated  with  respect  to  z  to  produce 
the  second  derivative.  The  terms 


d  u 

9  2 
z 


u 


and  — - 
9z^ 


kZ' 


kZ 


were  then  substituted  into  the  divergence  equa- 
tion. The  divergence  equation  then  became: 


(3) 


K 
kZ^ 


\  9x       9y  / 


Friction  velocity  was  obtained  by  integrating 
the  wind  profile  equation  for  fully  rough  surfaces. 
The  friction  velocities  are  then: 


|u|k 


{4a) 


In 


/Z  +  z  \ 


V .  = 


|v|k 


(4b) 


In 


/Z  +  z  \ 

(^) 


where  z^  is  the  surface  roughness  length  and  Z 
is  the  height  of  the  wind  above  ground. 

Rotational  flow  characteristics  were  intro- 
duced through  the  vorticity  equation.  Using  the 
same  assumptions  for  this  derivation  as  were 
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used  in  the  development  of  the  divergence  equa- 
tion, the  vorticity  equation  became: 


dv 


k'Vx 


dt 


8t 


9(0,  P*) 
'p  9(x,  y) 


(5) 


Again,  an  impulse  of  finite  time  interval  At^  was 
assumed,  so  that  the  mean  value  integration  gave: 


-(^+f)  5  - 


kZ 


(8v.       9u.\  r 
—  - - 
9x        8y  /         |_  ^ 


9(x,y) 


-  f6  - 


kZ 


(9v.  9u.\ 
—  ) 
9x        9y  /j 


The  divergence  and  vorticity  equations  are 
time-dependent  and  therefore  must  be  integrated. 
Divergence  was  integrated  directly  through  the 
mean  value  theorem  over  an  arbitrary  time  inter- 
val. Since  Little  is  known  about  the  rate  of  change 
of  divergence,  a  square  wave  impulse  of  interval 
At^  was  assumed.  The  integrating  time  factor 
then  became  YiAt^  when  the  magnitude  of  the 
acceleration  was  assumed  to  be  V2  the  observed  or 
maximum  value.  Divergence  was  then: 


[ 


9  =  hAtS      -c  0V  p.  -  c  V0-Vp, 


kZ 


\  9x        9y  /. 


(6) 


+  6 


Given  this  solution,  the  vorticity  equation  became 
the  ordinary  differential  equation  of: 


9t 


9(0, p..) 
+  7:6  =  -c    —  :  f6 


'P  9(x,y) 


K 
kZ^ 


\9x  "^97/ 


The  solution  is: 

-/6d 


C  =  e 


K 
kZ^ 


^9x  / 


9(0, P.,) 

■c  f6 

P  9(x,y) 


9y 


exp  (/6dt)dt 


-I-  ce 


-/6dt 


/I         X   ,      -5At  /2 
exp(iioAt  )  -f  ce 

The  arbitrary  time  intervals  for  the  divergence 
and  vorticity  were  defined  in  terms  of  phase 
velocities.  Because  a  rigid  lid  was  used  as  an  upper 
boundary  for  the  calculation  and  an  anelastic  ap- 
proximation was  used  for  compressibility,  the 
phase  velocity  for  divergence  was  taken  as  the 
speed  of  sound.  This  gave  a  time  interval  of: 

At  =  -g^^  ^ 
^6      300  m 

where  A  is  the  minimum  distance  between  compu- 
tational points.  Phase  velocities  for  vorticity  were 
taken  to  be  the  large-scale  windspeed,  so  that  the 
time  interval  for  vorticity  was: 


For  typical  distances  of  1  to  10  km  between  com- 
putational points  and  windspeeds  of  1  to  10  m  per 
second,  the  time  interval  is  of  the  order  of  1,000 
seconds.  Divergence  at  the  mesoscale  is  typically 
10"^ to  10"^ per  second.  These  magnitudes  sug- 
gested the  exponential  term  in  the  vorticity  equa- 
tion should  be  approximately  one,  except  under 
condition  of  light  winds  and  large  grid  spacing. 
The  constant  of  integration  was  then  taken  to  be 
the  large-scale  flow  vorticity.  The  vorticity  then 
became: 


9(0, P^) 


~c 


-  f6 


(7) 


.    P  9(x,y) 


kZ^      \  9x        9y  / 


Equations  (1)  and  (6)  for  the  frictionally, 
thermally  and  terrain-induced  flow  deviations  on 
divergence,  and  equation  (7)  for  the  frictionally 
and  thermally  induced  flow  deviations  on  rotation, 
were  then  used  to  calculate  the  wind  components. 


4 


Calculation  of  Local  Winds 

Wind  components  of  the  flow  were  calculated 
in  the  area  by  superimposing  the  terrain,  frictional, 
and  thermal  deviations  on  the  large-scale  wind. 
The  wind  deviations  were  obtained  by  separating 
the  flow  into  an  irrotational  divergent  component 
and  a  nondivergent  rotational  component.  The 
divergent  flow  was  described  with  a  velocity 
potential  and  the  rotational  component  with  a 
stream  function.  The  deviation  velocity  was  then: 

V  =  kxViJ;  + 

where  xp  is  the  stream  function  and  <|'  is  the  velocity 
potential.  The  vorticity  is: 


^  =  k'Vxv  =  k'VxkxVi|;  +  k'VxV({)  =  V  i|; 


and  the  divergence  is: 


(8) 


5  =  Vv  =  V'kxVip  +  V-Vc})  =  V  (1)  (9) 

Since  the  vorticity  and  divergence  were  known 
from  equations  (1),  (6),  and  (7),  the  stream  func- 
tion and  velocity  potential  could  be  obtained 
from  a  numerical  solution  of  the  elliptical  differ- 
ential equations  (8)  and  (9).  Decomposing  the 
velocity  into  the  west-east  component  and  the 
south-north  component,  and  superimposing  the 
background  flow,  the  velocities  are: 


and 


u  =  u    +  - 

o       dX  dy 


V  =  V     +  + 

o      dy  dx 


(10a) 


(10b) 


where  Uq  and  Vq  are  the  background  velocities. 


Flow  Near  the  Ground  Surface 

Airflow  along  the  terrain  surface  is  of  much 
more  interest  in  fire  or  air-quality  studies  than 
flow  at  some  fixed  elevation.  Since  the  equations 
require  calculations  to  be  made  on  a  flat  horizontal 
surface,  transformations  from  the  horizontal 
plane  to  the  terrain  surface  are  necessary. 

Data  required  by  the  model  are  observed  at  a 
fixed  height  above  the  ground  (generally  between 
1  and  3  m)  so  the  transforms  chosen  were  to  re- 
structure the  observed  data  to  a  flat  surface 
locally,  and  then  to  consider  the  model  as  being 


composed  of  a  series  of  adjacent  horizontal  sur- 
faces, each  of  shghtly  different  elevation  based 
on  the  elevation  of  the  local  terrain.  The  coordinate 
transform  required  expressing  the  derivations  of 
the  variables  observed  on  the  terrain  surface  as 
though  they  were  along  a  horizontal  surface.  If 
vertical  variations  of  the  variables  are  known,  the 
transform  for  any  quantity,  S,  is  then: 


9x , 


\      ^  ds\      _  IS  9h\ 

■/e      ^Vt  ~ 


where  the  unsatisfied  parentheses  indicate  deri- 
vations taken  along  either  constant  elevation  or 
along  the  terrain  surface. 

Second  derivations  were  obtained  similarly  as: 


=  -2 


9^  S 


9h 


9S  9^h 


9x  9z  9z      9z  di 


(ft)* 


9  2 
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The  higher  order  terms  involving  derivatives  of 
the  terrain  are  small  compared  to  the  Unear  term. 
Six  separate  derivatives  of  the  thermal  character- 
istics are  required  by  the  vorticity  and  divergence 
equations.  These  are  the  first  derivative  of  poten- 
tial temperature  and  pressure,  and  the  second 
derivative  of  pressure  in  the  cardinal  directions. 

Vertical  variations  of  these  thermodynamic 
variables  £U"e  required.  These  were  obtained  first 
for  potential  temperature  through  static  stability 
for  a  hydrostatically  balanced  atmosphere  as: 
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Vertical  decrease  of  pressure  was  represented 
by  hydrostatic  equilibrium: 
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Utilizing  these  relationships,  the  pressure 
derivatives  are: 
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Horizontal  derivatives  of  the  friction  velocities 
remain  unchanged,  since  the  model  is  designed 
for  use  in  the  constant  flux  layer  where  they  are 
constant  in  the  vertical.  This  shallow  layer  as- 
sumption then  gives: 
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These  linear  coordinate  transforms  neglect  the 
small  changes  in  the  grid  spacing  due  to  slope 
distances  projected  to  the  horizontal  grid. 


Data  Requirements  and 
Solution  Procedure 

Specific  data  required  by  the  model  are  tem- 
perature, pressure,  elevation,  and  roughness 
length  of  the  underlying  surface  for  all  computa- 
tional points.  Temperature  and  pressure  data  are 
transformed  thermodynamically  to  model  vari- 
ables by: 
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where  PqIs  a  reference  pressure  and  T  is  expressed 
in  an  absolute  temperature  scale.  While  any  con- 
sistent set  of  units  may  be  used,  the  appendix  list 
of  symbols  gives  constants  only  in  metric  units. 

Roughness  lengths  are  a  reasonably  constant 
value  describing  the  height  at  which  the  mean 
wind  flow  is  zero.  These  values  are  conveniently 
describable  to  acceptable  accuracy  by  classifying 
the  surface  into  such  categories  as  (but  not  limited 
to)  grass,  brush,  boulder  field,  and  so  forth.  These 
data  are  tabulated  in  numerous  texts  (Sutton 
1953,  Priestley  1959,  and  Sellers  1965);  a  sample 
Usting  is  provided  in  table  1  for  convenience  and 
completeness. 


Table  1 . --Re  1  at ionsh i p  between  surface  roughness  and  surface  features 


C 1  ass  i  f  icat  i  on  of 
underlying  surface 


Mud  flats,  ice 
Smooth  snow 
Snow,  pra  i  r 1 e 

Thin  grass  up  to  50  cm  high 

Thicl<  grass  up  to  50  cm  high 

Wheat 

City 

Orchard 

Fir  forest 


Typical  value  of 
roughness  length  (z.) 


0.001    (Sutton  1953) 
0.005  (Priestley  1959) 
0.01 

5  (Sutton  1953) 

9 

22  (Sellers  I965) 

165 
198 
283 
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Since  the  model  was  designed  for  use  in  mod- 
erately small  areas  of  finite  size,  several  single- 
valued  coefficients  are  required  to  complete  the 
model.  Most  important  of  these  are  the  static 
stabiUty  and  the  large-scale  wind  direction  and 
windspeed.  The  large-scale  winds  were  specified 
as  the  geostrophic  wind.  Less  important  coeffi- 
cients are  the  eddy  viscosity,  the  latitude  (if  the 
geographic  scale  of  a  particular  solution  is  large 
so  that  the  earth's  rotation  is  important),  and  the 
large-scale  vorticity  (if  strong  curvature  exists  in 
the  large-scale  flow). 

Wind  patterns  are  calculated  from  a  sequence 
of  approximations,  each  refining  the  previous  step 
by  introducing  a  more  refined  physical  process, 
and  superimposing  the  refinement  on  the  previous 
solution.  The  first  step  was  to  transfer  the  large- 
scale  background  wind  into  a  major  terrain- 
induced  modification  of  the  throughflow.  This 
step  provided  a  local  throughflow  wind. 

First,  the  divergence  was  calculated  through 
equation  (1)  for  a  smoothed  terrain.  This  smoothed 
terrain  was  obtained  by  averaging  the  four  sur- 
rounding points  with  each  computational  point; 
the  rigid  Ud  was  placed  well  above  the  terredn  in 
this  step.  Second,  the  terrain-induced  vorticity 
was  calculated  by  equation  (7)  with  the  thermal 
and  frictional  terms  neglected.  Once  these  vor- 
ticities  and  divergences  were  calculated,  they 
were  transformed  to  stream  functions  and  velocity 
potentials  through  numerical  solutions  of  equa- 
tions (8)  and  (9),  and  subsequently  resolved  into 
velocities  by  equation  (10).  These  first  approxi- 
mations to  the  wind  field  represented  the  influ- 
ences of  the  large  terrain  feature  on  the  through- 
flow. 

The  second  level  of  approximation  involved 
recalculating  the  terrain-induced  divergence  for  a 
shallow  layer.  Divergence  was  calculated  by 
equation  (1)  for  the  actual  terrain. 

The  next  step  required  using  the  large-scale 
terrain  flow  to  calculate  the  friction  velocities  by 
equations  (4a  and  4b)  and  making  the  coordinate 
transforms  of  equation  (11).  These  intermediate 
steps  provided  the  input  to  the  final  level  of  ap- 
proximation. 

Friction  terms,  defined  by  the  friction  velocity, 
and  the  spatial  derivations  of  temperature  and 
pressvire  after  the  coordinate  transform  were 
included  in  the  solution  of  equation  (6)  for  diver- 
gence and  of  equation  (7)  for  vorticity.  Background 
divergence  came  from  the  second  level  of  approxi- 
mation, and  the  background  vorticity  came  from 
the  divergence  and  vorticity  as  in  the  first  step 
through  numerical  solutions  of  equations  (8) 
and  (9).  Velocity  components  were  then  calculated 


through  equation  (10),  where  the  background 
velocities  were  defined  by  the  first  step. 

Evaluation  of  the  Model 

Validation 

There  are  several  philosophies  of  model  vali- 
dation, each  designed  to  accompUsh  specific 
objectives.  Because  this  model  is  intended  for 
use  in  applied  meteorology,  we  have  chosen  to 
use  actual  data  sets  from  mountain  mesoscale 
studies  to  vedidate  the  wind  field  estimates.  In 
adopting  this  approach,  we  have  imposed  a  severe 
test  since  the  model  is  limited  in  spatial  resolution 
by  computational  filtering  to  approximately  four 
times  the  grid  distance,  while  the  validation  data 
have  small-scale  spatial  resolutions  defined  by 
nearby  terrain  or  vegetation  and  sensor  response 
characteristics.  Also,  calculated  winds  represent 
an  average  wind  over  an  area  defined  by  the  grid 
mesh.  Such  a  test  is  desirable,  however,  because 
it  defines  operational  uncertainty  rather  than 
model  potential,  and  represents  actual  resolution 
attainable. 

Seven  data  sets  were  used  in  the  model  vali- 
dation. Six  came  from  studies  conducted  by 
Cramer  Mn  the  Oregon  Cascade  Mountains.  One 
data  set  came  from  central  Cahfornia  (Fosberg 
and  Schroeder  1966).  All  seven  depict  daytime 
conditions:  five  for  midafternoon,  and  two  for 
morning.  One  morning  case  occurred  under  over- 
cast skies  and  rain.  The  other  six  occurred  with 
strong  solar  heating.  Background  winds  ranged 
from  2  to  5  m  per  second.  The  Oregon  winds  were 
calculated  on  a  6  km  square  grid,  the  California 
case  on  a  16  km  square  grid.  Air  temperatures 
were  extrapolated  to  grid  points  from  the  field 
network  weather  stations. 

Evaluation  of  the  errors  in  both  windspeed 
(fig.  1)  and  wind  direction  (fig.  2)  demonstrates 
the  utihty  of  the  model.  The  root  mean  square 
windspeed  error  was  2.0  m  per  second.  Forty-five 
percent  of  the  calculated  windspeeds  fell  within 
1  m  per  second  of  the  observed  wind,  and  75  per- 
cent fell  within  2  m  per  second.  These  comparisons 
were  made  over  an  observed  range  of  0.5  to  10.5 
m  per  second.  Observed  wind  directions  were 
reported  on  a  16-point  compass  basis.  Therefore, 
the  error  analysis  is  presented  as  a  compass- 
point  error  rather  than  a  more  precise  eveduation 

^IVe  thank  Owen  P.  Cramer  of  the  U.S.  Department  of  Agri- 
culture, Forest  Service,  Pacific  Northwest  Forest  and  Range 
Experiment  Station,  Portland,  Oregon  for  unqualified  use  of 
original  data  for  the  model  validation  and  for  permission  to 
reproduce  unpublished  results  in  this  paper. 
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Figure  1. — Comparison  of  calculated  and 
observed  windspeeds.  Bold  diagonal  line 
defines  perfect  fit.  Forty-five  percent  of 
predictions  tall  within  1  m  sec~^  of  the 
observed  windspeed  and  75  percent  fall 
within  2  m  sec"^. 
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of  angle.  The  root  mean  square  direction  error 
was  1.9  points.  Sixty  percent  of  the  calculated 
directions  fell  within  1  compass  point  of  the  ob- 
served direction,  and  92  percent  fell  within  a 
2-point  envelope  of  the  observed  direction.  Ob- 
served directions  were  from  all  compass  points 
except  for  the  east  and  northeast  sections. 

Models  of  this  type  implicitly  require  a  num- 
ber of  specifiable  coefficients  which  tune  the  model 
to  particular  geographic  regions.  While  we  have 
attempted  to  specify  these  coefficients  through 
physical  arguments,  some  are  still  at  the  discre- 
tion of  the  user.  The  integration  coefficients  were 
totally  predetermined.  The  rigid  lid  coefficients 
for  potential  flow  were  selected  through  several 
criteria.  The  first  was  that  the  terrain  could  not 
penetrate  the  rigid  lid.  The  second  was  based  on 
a  large  number  of  solutions  where  the  Ud  was 
systematically  varied  until  a  high  level  of  sensi- 
tivity was  reached.  Lid  height  is  based  on  eleva- 
tion differences  across  the  grid  field.  The  user 
guideUnes  in  the  climatology  section  of  this  paper 
describe  the  techniques  for  determining  layer 
thickness.  They  are  cited  here  because  the  con- 
sequences are  reflected  in  the  error  analysis. 

Flow  Patterns 

Each  of  the  seven  cases  showed  similar  dy- 
namic processes  (table  2).  Thermally  and  terrain- 
produced  divergence  dominated  the  dynamics  of 
the  flow.  Divergence-generated  flow  perturbations 
were  2.5  to  4.3  times  as  large  as  vortically  gen- 
erated deviations.  Vortex  flow  accounted  for  11  to 
24  percent  of  total  kinetic  energy,  while  32  to 
69  percent  was  accounted  for  by  the  potential 
flow.  Strongest  local  wind  energetics  were  ob- 


served under  conditions  of  Ught  background  winds 
and  strong  surface  heating.  As  background  winds 
increased,  or  the  solar  heating  decreased,  the  local 
winds  became  a  weaker  component  and  less  impor- 
tant in  the  total  kinetic  energy  budget.  Tests  of 
background  windspeeds  for  several  hypothetical 
cases  beyond  the  range  of  data  used  in  validation 
show  that  at  0.5  m  per  second,  the  wind  patterns 
are  controlled  completely  by  the  thermal  patterns. 
At  10  m  per  second,  thermal  control  vanished  and 
terrain  influences  were  weak. 

Mesoscale  flow  characteristics  are  best  illus- 
trated through  a  detailed  examination  of  one  of 
the  case  studies.  The  area  is  in  the  Oregon  Cas- 
cades, south  of  Mt.  Hood,  north  of  Mt.  Jefferson, 
and  runs  from  the  Willamette  Valley  across  the 
crest  of  the  Cascade  Range  to  the  edge  of  the 
eastern  Oregon  plateau.  The  mountains  are  pene- 
trated by  a  major  drainage  on  the  west  side,  the 
Clackamas  River.  The  particular  case  selected, 
August  8,  1963  at  1600  PST,  was  chosen  as  rep- 
resenting a  typical  or  intermediate  situation.  The 
synoptic  weather  pattern  was  characterized  by  a 
weak  upper-level,  high-pressure  area,  and  west  to 
west-north-west  geostrophic  surface  flow  over  the 
Cascades.  Details  of  the  observation  and  analysis 
procedures  were  presented  in  Cramer  (1972)  and 
Cramer  and  Lynott  (1961,  1970).  Discussion  here 
is  limited  to  resultant  computations  and  estima- 
tion of  the  wind  field  from  the  analysis.  Only 
salient  features  of  Cramer's  analysis  are  presented 
here  for  completeness. 

Cramer's  data  analysis  of  temperature  and 
850  mb  surface  heights  were  interpolated  to  6-km 
grid  points.  These  data  were  converted  to  poten- 
tial temperature  and  normalized  station  pressures 
for  each  grid  point.  Terrain  elevation  at  each  grid 
point  was  interpolated  from  a  topographic  map. 


Table  2  . --Re  1  a t i onsh i p  between  flow  energetics  and  weather  features 


Date 

Loca  t  i  on 

T  i  me 

Flow  energy  accounted  for  by: 

Ratio  of 

Geostrophic 

Weather  features 

j  rrotat  iona 1 

w  i  nd  speed 

energy  to 

1 r  rota t  i ona 1 

Rotat 

ona  1 

rotat  iona 1 

LST 

energy 

energy 

energy 

percent 

m  per  second 

1963 

Aug . 

3  Oregon 
8 

0700 

60 

k 

Ik 

3 

2.5 

2.5 

Strong  morning  heating 

Aug . 

1600 

55 

1 

19 

8 

2.8 

2.5 

Strong  afternoon  heating 

Aug . 

19 

0700 

3'* 

2 

9 

1 

3.8 

k.3 

Weak  morning  heating-- 

extensive  clouds  and  rain 

Sept . 

k 

1600 

31 

5 

10 

7 

2.9 

h.S 

Strong  afternoon  heating 

Sept . 

5 

1600 

k2 

7 

15 

7 

2.7 

3.5 

Strong  afternoon  heating 

Sept . 

12 

1  300 

68 

5 

20 

7 

3.3 

2.5 

Strong  afternoon  heating 

1  962 

Strong  afternoon  heating 

Aug  . 

15  Calif. 

litOO 

k2 

3 

9 

.9 

'4.3 

3.0 
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The  underlying  surface  was  almost  exclusively 
characterized  by  conifer  forest,  so  the  roughness 
length  was  taken  as  a  constant  283  cm. 

Temperature  distribution  (fig.  3)  over  this 
terrain  generally  showed  a  sUght  decrease  with 
elevation.  West-facing  slopes  were  warmer  than 
east-facing  slopes.  The  major  river  valley  was 
much  warmer  than  surrounding  areas,  and  pro- 
vided strong  local  temperature  gradients.  In 
general,  the  area  east  of  the  Cascades  was  warmest 
as  a  result  of  morning  heating  and  topographic 
blocking  of  cool  marine  air  flowing  off  the  Pacific 
Ocean. 

Height  of  the  850  mb  pressure  surface  (fig.  4) 
generally  reflected  the  high-pressure  area  over 
northwest  Oregon.  Surface  heating  of  the  moun- 
tain slopes  produced  the  "thermal  low"  frequently 
observed  along  the  West  Coast.  This  "thermal 
low"  appeared  as  a  mesoscale  low-pressure  area 
through  the  higher  terrain  and  modified  the  large- 
scale  pressure  pattern. 

These  observed  data  were  transformed  to  the 
thermodynamic  variables  of  potential  temperature 
and  normaHzed  nondimensional  surface  pressure 
used  in  the  model  equations.  The  850  mb  surface 
height  was  reduced  to  the  terrain  surface  through 
the  hydrostatic  equation  using  the  static  stability 
from  the  Salem,  Oregon  rawinsonde  observation. 
Once  this  pressure  calculation  was  made,  the 
potential  temperatures  were  calculated  with  these 
pressures  and  the  interpolated  temperatures. 

Since  the  dynamic  influences  due  to  pressure 
gradients  are  completely  reflected  in  the  850  mb 
pressure  surface  height,  only  the  potential  tem- 
perature patterns  will  be  described  here  (fig.  5). 
Potential  temperature  showed  a  highly  organized 
pattern.  Higher  elevations  had  higher  potential 
temperatures.  Valleys  where  cooler  marine  air 
could  penetrate  from  the  coast  were  clearly  marked 
by  low  potential  temperature.  East-slope  high- 
lands had  highest  potential  temperatures. 

These  observed  fields  of  temperature  and 
pressure  were  used  to  calculate  the  fundamental 
kinematic  characteristics  of  divergence  and  vor- 
ticity.  Divergence  patterns  were  strongly  cor- 
related with  terrain.  Even  small  topographic 
features  such  as  narrow  ridges  and  valleys  stand 
out  in  the  divergence  patterns  (fig.  6).  Flow  diver- 
gence was  confined  to  lowlands  and  valley  bot- 
toms. Convergence  zones  were  on  ridgetops  and 
followed  the  terrain  contours.  Gentle  slopes  and 
large  continuous  upland  areas  on  the  east  slopes 
typically  had  weak  gradients.  Magnitudes  of  the 
divergence  were  typically  ±10"^  per  second.  The 
convergence  areas  in  regions  of  rugged  terrain 
were  much  smaller  than  divergence  areas,  and 
maximum  values  were  two  to  three  times  as  in- 


tense as  the  compensating  divergent  areas.  The 
intense  gradients  generated  in  the  rugged  ter- 
rain were  most  noticeable  between  the  Clackamas 
River  valley  and  the  surrounding  ridges.  Areas  of 
moderately  changing  terrain  had  divergent  and 
convergent  patterns  of  roughly  equal  area  and 
intensity. 

Vorticity  patterns  (fig.  7)  did  not  show  either 
the  response  to  small  terrain  features  or  the  large 
magnitude  found  in  the  divergence  fields.  Vor- 
ticity was  typically  ±10~5per  second  with  the 
strongest  areas  only  ±10~''per  second.  Ridgetops 
tended  to  have  anticyclonic  relative  vorticity  of 
10"'*  per  second  while  valley  bottoms  had  cyclonic 
vorticities  of  —  10"^per  second.  Areas  with  slowly 
changing  elevation  did  not  show  a  consistent 
terrain-vorticity  correlation.  In  general,  vorticity 
tended  to  have  a  much  longer  wave  length  than 
divergence.  The  only  generalizations  we  can  make 
about  the  patterns  are  that  major  ridges  tend  to 
favor  anticyclonic  vorticity,  while  lowlands  and 
major  valleys  tend  to  favor  cyclonic  vorticity. 
This  characteristic  was  probably  due  to  the  domi- 
nant influence  of  the  divergence  in  flow  properties. 

Resultant  wind  fields  (fig.  8)  showed  strong 
topographic  and  thermal  control  of  the  wind  pat- 
terns. Airflow  in  the  Clackamas  River  valley  and 
the  other  two  west-slope  watersheds  penetrating 
into  the  Cascades  showed  strong  topographic 
channeling  up  the  rivers  and  the  influence  of  slope 
winds  anghng  away  from  the  valley  axis.  West- 
slope  areas  away  from  these  major  drainages 
generally  showed  upslope  flow,  even  on  the  lee 
sides.  Strongest  winds  were  calculated  on  high 
ridgetops  and  at  the  heads  of  canyons  as  they 
reached  the  highlands.  Lee  effects  were  most 
noticeable  along  the  eastern  slope.  Winds  flowing 
across  the  ridge  in  the  northeast  portion  of  the 
area  were  reduced  from  5  m  per  second  to  1  m  per 
second  by  the  upslope  influences  on  the  lee  side. 
Strong,  thermally  driven  winds  were  calculated 
in  the  south-central  and  southeastern  portion  of 
the  area,  where  prevaihng  winds  were  reduced  and 
influences  of  the  penetrating  west-slope  canyons 
were  negligible.  Those  variables  which  had  great- 
est impact  on  model  reUability  were  determined 
by  sensitivity  analysis.  Both  processes  and  vari- 
ables were  systematically  excluded  and  error 
analysis  conducted  on  all  seven  cases.  This  analy- 
sis indicated  that  frictional  influences  and  the 
vorticity  equation  affected  the  resultant  wind 
field  only  slightly.  Terrain  and  thermal  influences 
totally  dominated  the  calculations.  Neglect  of 
either  created  large  errors  in  the  wind  field  in  the 
test  cases.  Further  analysis  of  these  processes 
indicated  that  pressure  variation  could  be  ex- 
cluded with  introduction  of  a  small  error,  but  that 
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Figure  3. — Temperature  field  for  1600  PST,  August  8,  1963.  Temperatures  are  m  degrees  k;  warmest  areas  are  on  the  east  slopes 

of  the  Cascade  Mountains. 


Figure  4. — Height  of  850  mb  surface  for  1600  PST,  August  8,  1963.  Height  is  in  meters  above  m.s.l.  Thermal  trough  is  oriented 

southwest  to  northeast- 


Figure  5. — Distribution  of  potential  temperature  for  1600  PST,  August  8,  1963.  Temperatures  are  In  degrees  k.  Inversion  of  marine 

air  into  the  mountains  Is  confined  to  the  Clackamas  River  valley. 
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Figures. — Horizontal  divergence  field  tor  1600  PST,  August  8,  1963.  (Units  are  10  ^sec~i.)  Strongest  gradients  of  divergence  are 

between  the  west  side  drainages  and  the  surrounding  ridges. 


Figure  7.— Vertical  component  of  vorticity  for  1600  PST,  August  8,  1963.  (Units  are  10'^sec  M  Anticyclonic  vorticity  is  on  ridgetops 

and  cyclonic  vorticity  is  in  valleys. 


Figure  8. — Calculated  and  observed  wind  field  for  1600  PST,  August  8,  1963.  Arrows  indicate  direction  of  wind.  Speeds  are  given 
by  sum  of  barbs  and  pennants  on  arrow  tails.  Pennants  indicate  5  m  sec-i;  a  full  barb,  2  m  sec"i;  a  half  barb,  1  m  sec~i.  Calculated 
wind  field  is  in  red;  observed  wind  field  is  in  black. 
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the  topographic  and  temperature  field  must  be 
represented  in  complete  detail  for  the  model  to 
provide  reaUstic  and  verifiable  results. 


Topography,  Climatology,  and  Winds 

Two  parameters  are  required  if  the  model  is 
to  define  the  topographic  influence  on  through- 
flow.  The  first  of  these,  the  depth  of  the  rigid  lid 
at  the  top  of  the  slab,  should  normally  be  expected 
to  be  1,500  to  2,000  m  above  the  spatially  smoothed 
terrain.  The  second  parameter,  the  thickness  of 
the  shallow  layer,  refines  the  first  estimate  of 
terrain  influences.  It  is  1.25  times  the  difference 
between  the  smoothed  terrain  and  the  actual  ter- 
rain. These  guidelines  for  the  relationship  between 
throughflow  and  terrain  were  determined  from  a 
large  number  of  test  cases  where  the  coefficients 
were  varied  systematically.  The  final  choice  was 
a  judgment  of  which  relationships  produced  the 
most  acceptable  resultant  wind  field. 

Grid  point  spacing  is  also  related  to  topo- 
graphic variations.  In  rugged  terrain,  a  closely 
spaced  grid  is  required  to  resolve  the  local  winds, 
while  in  gently  rolling  terrain,  a  low-resolution 
grid  will  give  acceptable  wind  fields.  Guidelines 
for  selecting  appropriate  grid  spacing  were  de- 
veloped from  procedures  similar  to  those  used  in 
evaluating  layer  thicknesses. 

Numerous  analyses  were  carried  out  with  the 
same  data  sets  so  that  performance  could  be 
evaluated  by  overall  resultant  wind  field  errors. 
Grid  spacings  between  1  and  6  km  were  evaluated 
for  maximum  terrain  differences  of  300  to  1,500  m. 
These  analyses  indicated  that  grid  spacing  should 
be  reduced  by  1  km  for  every  increase  in  terrain 
range  of  300  m.  If  the  terrain  within  an  area 
10,000  km  square  varies  by  300  m  or  less  between 
the  highest  and  lowest  spatially  averaged  eleva- 
tions, then  a  6-km  grid  would  be  appropriate.  If, 
on  the  other  hand,  the  terrain  variation  were 
1,500  m  a  1-km  grid  would  be  appropriate.  These 
spacings  must  be  modified  to  the  smallest  grid 
spacing  if  the  wind  field  is  intended  to  be  used  for 
fine-scale  analysis  of  pollution  patterns  or  fire 
spread,  or  to  larger  spacing  if  broad  geographic 
coverage  is  required.  Grid  spacing  of  500  m  is 
probably  the  minimum  the  model  will  treat  prop- 
erly. Close  grid  spacing  does  not  permit  a  50-  by 
50-point  grid  to  include  major  terrain  features. 
These  two  nonmodel  dependent  factors  (terrain 
complexity  and  computer  resources)  must  be  con- 
sidered in  selecting  the  grid  interved.  The  spacings 
specified  here  represent  maximum  spacing— not 
optimum  for  any  particular  appUcation. 


The  seven  cases  in  this  analysis  all  came 
from  detailed  mesometeorological  research  studies. 
Such  intensive  data  would  not  normally  be  avail- 
able. Detailed  temperature  patterns  can  be  esti- 
mated from  chmatological  relationships  developed 
for  slope  temperatures.  Super  and  Grainger  (1969) 
and  Hayes  (1941)  made  detailed  studies  of  slope 
temperatures  by  aspect  and  time  of  day.  Baker 
(1944),  using  maximum  and  minimum  tempera- 
tures, characterized  the  slope  temperatures  for 
20  mountainous  areas  of  the  Pacific  and  Rocky 
Mountain  States.  These  chmatological  relation- 
ships will  approximate  the  temperature  distribu- 
tion when  the  available  observations  and  topog- 
raphy are  used  to  fine  tune  the  slope  functions. 

Using  the  Model  for  Planning — A 
Case  Study 

Consider  a  case  where  engineers  and  architects 
are  designing  a  new  plant  at  a  mine  mouth.  The 
mineral  extraction  process  will  create  a  large 
volume  of  tailings  over  a  period  of  years.  These 
tailings  could  be  hot.  In  the  design  study  the 
questions  are  asked:  (1)  What  will  be  the  effect  of 
changing  the  terrain  in  the  immediate  area  by 
filUng  the  gulches  and  small  valleys  with  tailings 
on  the  downwind  dispersal  of  stack  contaminants?, 
and  (2)  What  will  be  the  effect  on  downwind  dis- 
persal if  the  tailings  are  approximately  70  °C 
warmer  than  the  surrounding  areas? 

The  wind  model  may  be  used  to  estimate 
probable  consequences  of  the  mining  activity  on 
pollution  dispersal  patterns.  To  illustrate  this 
use,  we  assumed  a  Ught  wind  of  0.45  m  per  second 
(1  mile  per  hour)  from  the  southwest;  temperature 
distribution  was  given  by  climatology  for  night- 
time conditions.  Wind  patterns  before  construc- 
tion showed  well-organized  drainage  flow  into 
and  down  the  valleys  and  gulches  (fig.  9). 

Two  cases  of  land  change  were  simulated  to 
anticipate  influence  on  the  local  wind  patterns.  In 
the  first  case,  the  valleys  and  gulches  were  filled 
with  heated  tailings.  Changes  in  the  flow  were 
calculated  near  the  modified  landform.  Heated 
areas  increased  the  windspeeds  and  changed  the 
directions  to  a  uniform  south-southwest  flow 
down  the  major  gulch  (fig.  10).  Airflow  in  the 
large  river  drainage  to  the  north  was  turned  and 
flowed  up  canyon.  The  second  case  considered 
was  that  the  same  gullies  and  canyons  were  filled 
with  tailings,  but  the  taiUngs  were  not  heated. 
This  modification  also  increased  and  organized 
the  drainage  flow  (fig.  11).  Winds  were  from  the 
southwest  over  the  tailings  for  this  case.  How- 
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indicates  plant  and  mine  site. 


Figure  1 0. — Simulated  wind  field  for  a  case  with 
heated  tailings. 


Figure  9. — Local  wind  pattern  before  construction 
of  a  hypothetical  plant  and  mine  tailings. 


indicates  plant  and  mine  site. 


Figure  1 1 . —  Simulated  wind  field  for  a  case  with 
cold  tailings. 


indicates  plant  and  mine  site. 
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ever,  air  did  not  flow  up  the  canyon  from  the  major 
river  drainage.  The  two  cases  predicted  increases 
in  windspeed  of  as  much  as  1  m  per  second  and 
direction  changes  of  180  degrees  near  the  planned 
mine  activities. 

Summary  and  Conclusions 

Local  thermally  driven  flows  such  as  moun- 
tain and  valley  winds  are  simulated  reasonably 
well.  Interactions  of  these  local  flows,  and  influ- 
ence of  the  larger  scale  potential  flow  on  these 
local  winds,  are  also  included  in  the  resultant 
flow.  These  thermally  driven  flows  Eire  determined 
properly  primarily  because  the  thermal  forcing 
processes  are  stationary  and  locked  to  the  terrain 
features.  Also,  these  thermal  features  are  in  a 
quasi-steady  state  for  several  hours.  Thus,  neglect 
of  the  advective  process  is  permissible,  and  the 
impulse  integration  to  steady  state  approximates 
the  flow  field. 

Thermal  modification  of  potential  flow  across 
the  terrain  features  dominates  the  driving  forces. 
Potential  flow  was  three  times  as  important  as 
rotational  flow. 

Error  analysis  of  seven  case  studies  and  the 
apphcation  example  indicated  that  the  model 
could  be  used  as  an  aid  in  fire  and  air-quality  pro- 
grams because  of  the  low  data  requirements  and 
the  reduced  computer  effort  required. 

The  wind  model  developed  here  was  based  on 
simplifications  of  the  complete  set  of  equations 
governing  atmospheric  flow  over  complex  terrain. 
These  simpUfications  imposed  several  Umitations 
on  the  use  of  the  model,  both  in  types  of  flow  that 
can  be  depicted  and  in  spatial  coverage  of  the 
solutions.  Moving  flow  systems  and  flows  which 
are  characterized  by  intermittency  or  by  separa- 
tion are  explicitly  excluded,  because  (1)  advection 
was  eUminated,  and  (2)  the  model  is  intended  to 
depict  quasi-steady  state  processes.  Coefficients 
describing  static  stabiHty  and  overall  wind  through- 
flow  place  spatial  limits  on  the  model  because 
(1)  constant  static  stability  prohibits  larger  scede 
baroclinicity,  and  (2)  the  single-valued  wind  vector 
prohibits  large-scale  flow  discontinuity.  Also,  the 
neglect  of  advective  processes  excludes  strong 
interaction  of  the  mesoscale  with  larger  scales, 
and  does  not  allow  for  modification  of  the  down- 
wind enviroiunent. 
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APPENDIX 
List  of  Symbols 


Symbol 


Definition 


Units 


f 

g 
H 
h 
k 

k 
K 
P 
Po 


R 
T 

t  _ 

U,  U,  Uq 
V 

V,  V,  Vo 
V* 

X 

y 

z 

0 

p 

a 

V 

A 

At^ 


Specific  heat  at  constant  pressure 

Coriolis  parameter 

Acceleration  of  gravity 

Thickness  of  rigid  lid  layer 

Thickness  of  second  approximation  to  flow 

von  Karman  constant 

Unit  vector  in  vertical  direction 

Turbulent  viscosity 
Pressure 

Reference  pressure 
Non-dimensional  pressure 


R 


Ideal  gas  content 
Temperature 
Time 

Windspeed  in  west  to  east  direction 
Friction  velocity  in  west-east  direction 
Horizontal  velocity  vector 
Windspeed  in  south-north  direction 
Friction  velocity  in  south-north  direction 
Vertical  windspeed 

Horizontal  coordinate  in  west-east  direction 
Horizontal  coordinate  in  south-north  direction 
Vertical  coordinate 
Surface  roughness  length 
Horizontal  divergence 
Vertical  component  of  vorticity 

T 

Potential  temperature  ^  = 


Density  of  air 
Static  stability  a 


1 


30 
dz 


Velocity  potential 
Stream  function 

Horizontal  vector  gradient  operator 
Horizontal  laplacian  operator 
Horizontal  grid  spacing 
Impulse  time  for  vorticity 
Impulse  time  for  divergence 


1.003  joules  gm 

sec-i 

9.8  m  sec~2 

m 

m 

0.4 

dimensionless 

m^  sec~i 
dynes  cm -2 
10^  dynes  cm" 


-1  o 


K" 


0.287  joules  gm  '  °K" 
DegK 
Sec 
m  sec 
m  sec 
m  sec 
m  sec 
m  sec 
m  sec 
m 
m 
m 
m 

sec-i 
sec-i 

DegK 

kg  m~3 

m^  sec 7^ 
m^  sec~i 
m-i 
m-2 

m 
sec 
sec 
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